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We compare the physical potential Vd(R) of an external quark-antiquark pair in 

c3 



the representation D of SU(N), to the color-Coulomb potential V CO ui(R) which is the 
instantaneous part of the 44-component of the gluon propagator in Coulomb gauge 
^| Da{x,t) = V cou \(\x\)5(t) + (non-instantaneous). We show that if Vd(R) is confining, 



X 



liniR^oo Vd(1?) = +oo, then the inequality Vd(R) < —CdV cou \(R) holds asymptotically 
at large R, where Cd > is the Casimir in the representation D. This implies that 
—Vcoui(R) is also confining. 



1. Introduction 



The problem of confinement of color charge has been with us for a long time. There 
are many approaches to this problem such as dual Meissner effect by monopole condensa- 
tion [1], effective string theory [2], center dominance [3], and color-Coulomb potential [4]. 
One seeks to choose variables so that the most important degrees of freedom have a simple 
expression. For this purpose gauge fixing can be a useful technique. 

Confinement is most commonly characterized by the behavior of Vd(R), the gauge- 
invariant potential energy between an external quark pair at separation R in the rep- 
resentation D of the gauge structure group SU(N). It may found from a rectangular 
Wilson loop Wd(R, T) in the representation D, but for our purposes it is more conve- 
nient to obtain it from the correlator, (P£>(x)P^(y)) , of a pair of Polyakov or thermal 
Wilson loops at x and y in the representation D, on a Euclidean lattice of period T 
in the 4-direction. The Polyakov loop is the lattice analog of the continuum expression 
Pd(x) = tr[Pexp(J Q i Bi4 (x, t)dt], where A Dyfl = A a ^t a D , and the t a D satisfy the Lie al- 
grebra commutation relations \t a Dl t h ^\ = f abc t c D in the representation D. As discussed 
recently [2], this correlator has the expansion 



where the E n ^^ are the eigenvalues, H^> n = E n ^^ n , of the lattice QCD hamiltonian 
if, specified below, that includes an external quark and anti-quark at x and y in the 
representation D. In the large-T limit, the sum is dominated by the first term, with lowest 
energy eigenvalue, Eq^^. It is rotationally symmetric in the continuum limit, and we 
identify the physical quark-antiquark potential Vd(-R) with this energy eigenvalue, after 
separation of divergences, 



Here Eq is the energy of the vacuum state in the absence of an external quark pair, and A 
is the diverging self-energy of the external quarks. According to the Wilson confinement 
criterion, which is expected to hold in pure gluodynamics without dynamical quarks, but 
with external quarks in the fundamental representation, Vp(R) diverges linearly at large 
.R, Vp (R) ~ crR where a is the conventional string tension. 




(1.1) 



E 0i \ Sl _ Sa \ =E + A + V D (\x-y\). 



(1.2) 
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In the Coulomb gauge, there is a simple scenario [4] that attributes confinement of 
color charge to the long range of the the color-Coulomb potential, V cou \(R). This quan- 
tity characterizes the instantaneous part of the 44-component of the gluon propagator 
(Al(x,t)A b 4 (0,0)) = D M (x,t)5 ab = ^coui(N) $(t) S ab + (non-instantaneous). Since A 4 
couples universally to the color-charge, this can account for confinement of color-charge, 
provided that V CO ui(R) is indeed long range. A remarkable feature of the Coulomb gauge 
in QCD, a property not shared by any Lorentz gauge, is that A4 = goA^ = g r A^ 
is a renormalization-group invariant [4], [5]. Here go and A± , and g r and A£ are, 
respectively, the unrenormalized and renormalized charges and perturbative gauge con- 
nections. This means that -D44, and hence also its instantaneous part V^ ou i(i2), is in- 
dependent of both the cut-off A and the renormalization mass (i. This property allows 
the fundamental QCD quantity, V cou i(R), the instantaneous part of the gluon propaga- 
tor, to be identified with the phenomenological QCD potential [6] and [7]. Its fourier 
transform, KouiG^I), provides a convenient definition of the running coupling constant, 
a s (|£|/A coul ) = S c 2 ou1 (|£|)/(4tt) = k 2 V coul (\k\)/(4TTX ), where x = 12N/(11N - N f ), A coul 
is a finite QCD mass scale, and Nf is the number of quark flavors [8]. The result obtained 
here means that if the Wilson criterion for confinement is satisfied, then, with this defini- 
tion, the running coupling constant a s (|/c|/A cou i) diverges in the infrared like 1/fc 2 , a clear 
manifestation of infrared slavery. 

We shall show that a necessary condition for confinement according to the Wilson 
criterion is that the instantaneous color-Coulomb potential be confining. In symbols: if 
liniR^oo Vd (R) = +00, then Vd(R) < — Co V CO ui(R) holds asymptotically at large R. 
The minus sign occurs because antiquark has opposite charge to quark. Here Co = 
— Ylai^h) 2 > is the (positive) value of the Casimir invariant, and in the fundamental 
representation Cp = (N 2 — 1)/(2N). A considerable simplification is hereby achieved 
because Vd(R) is defined by means of a path-ordered exponential that involves gluon n- 
point functions of all orders, whereas V cou i(R) is defined in terms of the gluon 2-point 
function, -D44. We also note the striking numerical result [9], [10] that Vo(R) exhibits 
Casimir scaling, Vd(R) = (Cd/Cf)Vf(R) quite accurately at least in a rather large range 
of r and 8 representations D. This suggests that the above bound may be saturated in this 
range, for this would explain Casimir scaling, that is not easy to understand otherwise [11]. 
The result also makes it imperative to extend present programs to calculate V CO ui(R) 
numerically [12], and analytically from first principles [7], and to derive phenomenological 
quantities from it [6]. 
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2. Lattice Coulomb-gauge QCD hamiltonian 

The energy £?o,|5i-x 2 | ^ s °f course g au g e invariant, and the lattice QCD hamiltonian 
H may be chosen in any gauge. Its most familiar form is in the temporal gauge U4 = 1, 
corresponding to A 4 = 0, 

tftemp = gl&a)- 1 + 2(g 2 a)- 1 J2^ ^U p , (21) 

x,i p 

where J2 P * s the sum over all spatial plaquettes p (on a single time-slice). Here £^,i is 
the color-electric field operator that satisfies [^g^U^j] = it a Ujjj8g,ySij and [£^ v £^j] = 
—if ahc 8x,y^ij£$i- We place an external quark at x\ in the representation D, and an 
external anti-quark at xi in the representation D* , with color vectors that act on the first 
and second indices of the wave-functional ^ a p(U) according to (A"\I/) Q/ 3 = (A^,) a7 \l/ 7i a and 
(A2^/) a /3 = — (^£>)/? 7 *a7) where \ a D = it a D . In the temporal gauge, the color charges of 
the external quarks do not appear in the hamiltonian H temp , but rather in the subsidiary 
condition G a (x)^ = 0. This is an expression of Gauss's law, for G a (x) is a precise lattice 
analog of the continuum Gauss's law operator G a (x) = — (D-E) a (x) + p^ u (x), where 
Ef(x) = i 5A i^ , Df c = 5 ac di + f ahc A\ is the gauge-covariant derivative, and p^ u (x) = 
\i5(x — xi) + \25(x — X2) is the color-charge density of the external quarks. In the temporal 
gauge G a (x) is the generator of 3-dimensionally local gauge transformations of the quark 
and gluon variables, satisfying [G a (x), G b (y)] = i5(x — y)f abc G c {x), and the subsidiary 
condition is the statement of gauge invariance of the wave functional. 

One would expect that Gauss's law is essential for confinement, and the lattice 
Coulomb hamiltonian H coxl \ [13] may be derived from H temp by solving Gauss's law as 
subsidiary condition [14] . For our purposes the resulting lattice Coulomb hamiltonian has 
the same structure to the continuum Coulomb hamiltonian [15]. To simplify the exposi- 
tion, we shall use continuum language, but it is understood that this is short-hand for the 
correct lattice kinematics, and divergences are controlled by use of the lattice Coulomb 
hamiltonian, as will be made clear. 

To get to the Coulomb gauge from the temporal gauge, one integrates out the gauge 
degrees of freedom using the Faddeev- Popov formula in all gauge-invariant matrix elements. 
In particular for the hamiltonian, one obtains H cou \, defined by its matrix elements 

(*i,ifcoui* 2 ) = ^cL4 tr detM (1/2) J d 3 x \ S *(E?* 1 )*E?* 2 + g^B 2 **], (2.2) 
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there the wave-functionals ty a p(A tr ) depend only on 3-dimensionally transverse continuum 
configurations diAf = 0, and a contraction on color indices is understood. The color- 
magnetic field is given by Bf = d 2 Af a - d s Af a + f ahc Af b Af c , etc., and the color- 
electric field by Ef = Ef ,a — di<j) a , where Ef ,a = i t-^e , and <p a (x) is the color-Coulomb 

i 

potential operator. In this matrix element, <f) a (x) acts directly on the wave functional 
The definition of H cou \ is completed by specifying that cf) a (x)^ = (M ~ x p p h ys ) a (x)^ , which 
expresses (p a (x)fy in terms p qu and transverse gluon variables only. This is the solution 
of the subsidiary condition G a (x)^ = 0, or M ac (A tr )(j) c ^ = Pp hys ^. Here M ac (A tr ) = 
-Df c (A tT )di = -d t D? c (A tT ) = -d 2 5 ac - f abc Af' b d l is the 3-dimensional Faddeev-Popov 
operator, and Pp hys = —f ahc A^'Ef ,c + p® u is the color-charge density of the external 
quarks plus the color-charge of the dynamical gluon degrees of freedom only. The associated 
color charge, Q a = J d s x Pp hys (x) may be identified with the physical color charge, for it 
generates global gauge transformations on all variables [Q a 7 Af' b ] = if abc Af' c etc., and 
satisfies [Q a , Q b ] = if abc Q c . The second term in M is characteristic of non-Abelian gauge 
theory. It is responsible for anti-screening because, for typical configurations, this term 
produces a small denominator in M~ x . The subscript A on the integral j A cL4 tr means 
that a region that includes only one Gribov copy is integrated over. This may be chosen 
as in the minimal Coulomb gauge, but the proof does not depend on the particular way 
this is chosen. 

3. Bound on Vr>(R) from trial wave function 

The energy E\ Sl _^ 2 \ = (\&, H cou \^>) of any trial wave function ^> provides an up- 
per bound on the ground-state energy, E Q ^ 1 _^ 2 \ < E^_^ 2 \. As trial function we take 
the product wave function, # Q/3 (A tr ) = N D 1/2 5 af3 $ (A tr ). Here $ (A tr ) is the exact 
wave functional of the vacuum state in the absence of external quarks, and A"^ 1 ^ 2 ^^, 
where Njj is the dimension of representation D, is the external quark-pair state of to- 
tal color-charge zero, (Ai + \2) a ^ = 0. The Coulomb hamiltonian has the decomposi- 
tion H cou \ = Hg\ + H g i qu + H qu qu , that follows from the decomposition of the color- 
charge density p^ hys = -f<*c£*> b Ef' c + p a qu in * = M _:l pphys Here H gi is the 
Coulomb hamiltonian in the absence of external quarks, i?gi,q U is linear in p qu , and 
-Hqu,qu = (1/2) / d 3 x (<9iM -1 p qu ) 2 (if). In the last expression there is no ordering prob- 
lem because p qu commutes with A tr . By definition of $o we have H g \$o = E &o, where 
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Eq is the vacuum energy in the absence of external quarks. From (\1>, if g i\l/) = E , and 
Hgi j(lu ^) = 0, we get for the trial energy, 

E Wl _ 3al = E + A' - C D (d> , [M-'i-d^M-^J^o) (3.1) 

(no sum on a), where we have used (AJA£) = (A^ 2 - l)" 1 5 ab (AJ A§) = — (A^ 2 - 
l^d^Co. Here A' is another self energy of the external quarks that is indepen- 
dent of x\ and X2- The inequality E ^ 1 _^ 2 \ < E\$ 1 _£ 2 \ reads A(A) + Vd(R, A) < 
A'(A) - Cd($o, [M ~ 1 (-9 2 )M~ 1 ]^ i ^ 2 $ )- We have cancelled the vacuum energy E 
that diverges with the volume of space and, having done so, we may take the volume of 
space to infinity, keeping the ultraviolet cut-off A = a -1 in place, where a is the lattice 
spacing. Here A(A) and A'(A) are self-energies of the external quarks. The formula [8] for 
the color-Coulomb potential, Kouid^i - x 2 \)d ah = ($ , [M-^-d^M" 1 ]^ ^ $ ), allows 
us to write the inequality as A (A) + V D (R, A) < A' (A) - C D V cou i(R, A). 

We have inserted a dependence on the cut-off A in Vd(R, A) and (i?, A) , be- 
cause these are lattice quantities that depend on the lattice spacing, a = A -1 . How- 
ever, having separated out the self-energies, both the quark potential Vd(R, A) and the 
lattice color-Coulomb potential V cou \(R, A) have finite, A-independent continuum limits, 
li m A-><x> Vd(R, A) = V D (R) and lim A ^oo Koui(-R> A) = V cou i(R). For V D (R) is a physi- 
cal energy, and V cou \(R) is a renormalization-group invariant, as noted in the Introduc- 
tion. If Vd(R, A) is confining, lim^^oo Vr> (R, A) = +oo, then, for sufficiently large i?, 
the self-energies A(A) and A'(A) are negligible compared to Vd(R, A), and the inequality 
Vd{R, A) < — (7d V^ouK-R, A) holds asymptotically at large R, for finite cut-off A. This 
bound also holds in the continuum limit, because dimensional and renormalization-group 
considerations tell us that the terms that vanish as A — > oo are of relative order l/(Ai?) n , 
where n is positive, so they also vanish asymptotically at large R. 1 We conclude that, if 
Vd(R) is confining, lim ji^ oq Vd(R) = +oo, then in the continuum limit, the inequality, 
Vd(R) < — Cd V cou \(R), holds asymptotically at large R, as asserted. 

1 This condition is necessary, as shown by the following counter-example. Take Vd(R, A) = aR, 
and y coul (i?,A) = c/R + m i R 2 /A, with self-energies A(A) = a A, and A'(A) = (a + 1)A. The 
inequality at finite A reads aR < A + c/R + m 4 R 2 /A = 2m 2 R + c/R + (A - m 2 R) 2 /A. It is 
satisfied for all finite A and R, provided that a < 2m 2 and c > 0. But the continuum limit of 
V cou \{R) is c/R, and aR < c/R is not satisfied at large R. 
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4. Conclusion 

The bound implies that if the potential between external quarks in the fundamental 
representation increases linearly at large R, Vp(r) ~ aR, where a is the standard string 
tension, then the color-Coulomb potential V cou i(R) increases at least linearly at large R, 
and moreover if its increase is also linear, — V cou \(R) ~ cr cou i-R, as has been conjectured 
[4], where a cou \ is a string tension that characterizes V cou i(R), then the conventional string 
tension satisfies the bound a < (N 2 — 1) / '(2N)a cou \. 

What has been learned about QCD dynamics? We have found that if the Wilson 
confinement criterion holds, then V cou i(R), the instantaneous part of the gluon propagator 
D44 in Coulomb gauge, is confining. Moreover, from V cou \(\x — y\) = ([M _1 ( — d 2 )M~ 1 ]^^), 
this can happen only if the Faddeev- Popov or ghost Green function, [M~ 1 (A tr )]^^ : is long- 
range for configurations A tr that dominate the functional integral. This confirms the con- 
finement scenario originally proposed by Gribov [16], and advocated by the author [4]. [The 
scenario reads, in brief, that in the minimal Coulomb gauge, configurations are restricted 
to the Gribov region, where the Faddeev-Popov operator is positive, M(A tr ) > 0. The 
boundary of this region occurs where the lowest eigenvalue of M(A tr ) vanishes, \o(A tr ) = 0. 
Moreover the dimension n of configuration space is very large, being of the order of the 
volume V of the lattice. Entropy favors a population highly concentrated close to this 
boundary, where \o(A tr ) is small, for the same reason that, in a space of very high dimen- 
sion n, the density of a ball r < ro is very sharply peaked near ro, being given by r n ~ 1 dr. 
Consequently, for the configurations that dominate the functional integral, M~ 1 (A tr ) is 
enhanced, and thus also V cou \(R).] 

To simplify the exposition, we considered gluodynamics without dynamical quarks, 
but the proof holds if they are included. However if dynamical quarks are present in the 
fundamental representation F, as occurs in nature, then the physical potential Vp(R) be- 
tween external quarks does not manifest confinement. For at some radius Rb the string 
breaks by polarization of sea quarks from the vacuum, and for R > R^, Vp{R) represents 
a residual potential between a pair of mesons, analogous to the van der Waals potential. 
In this case the bound obtained here does not imply that V COVl \(R) is confining. Neverthe- 
less, according to the confinement scenario in Coulomb gauge, V CO u\(R) is a fundamental 
quantity that remains linearly rising even when Vf(R) is not. It is precisely the linear rise 
of V cou \(R) that causes string- breaking, by making it energetically preferable to polarize 
sea-quarks from the vacuum. 
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